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ABSTRACT 



We employ N-body, smoothed particle hydrodynamical simulations, including detailed treat- 
ment of chemical enrichment, to follow a gas-rich merger which results in a galaxy with disk 
morphology. We trace the kinematic, structural, and chemical properties of stars formed be- 
fore, during, and after the merger. We show that such a merger produces two exponential disk 
components, with the older, hotter component having a scale-length 20% larger than the later- 
forming, cold disk. Rapid star formation during the merger quickly enriches the protogalactic gas 
reservoir, resulting in high metallicities of the forming stars. These stars form from gas largely 
polluted by Type II supernovae, which form rapidly in the merger-induced starburst. After the 
merger, a thin disk forms from gas which has had time to be polluted by Type la supernovae. 
Abundance trends are plotted, and we examine the proposal that increased star formation during 
gas-rich mergers may explain the high a-to-iron abundance ratios which exist in the relatively 
high-metallicity thick disk component of the Milky Way. 
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Introduction 



That major mergers of two gas-rich disk galaxies can sporn a disk galaxy was shown in lSpringel fc Hernguist 

(j2005l ). Disk galaxies can result from such events for a large range of mass ratios, orbits, and rotational 



velocities of progenitor galaxies ( Robertson et al. 2006[). The discovery that thick disks are common, and 
perhaps even ubiquitous in disk galaxies |Dalcanton fc BernsteinlboodlMouldboOsI ). and their old age means 
that the thick disk holds important clues to the formation of such galaxies. It also implies t hat disk galaxies 



formin g from a major merger should have a thick disk component. Examining Figure 2 of [Robertson et al 



(|2006[ ). it is clear that a relatively hot (thick) disk component will naturally form in such mergers, motivat- 
ing our desire to further investigate such events. In general terms, mergers were presumably increasingly 
gas-rich at increasingly early epochs. The possibility th at gas-rich mergers play a n essential role in disk 
galaxy formation has gained momentum in recent years. iBrook et all (|2004bl liooil suggest that the thick 
disk is formed in a high-redshift, gas-rich merger epoch, characteristic of a ACDM Universe, during which 
the galaxy is born. High-redshift mergers being gas-rich is also consistent with other constraints including a 
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low-mass, metal-poor stellar halo (Brook et al.ll2004al) . the chemical abundanc e patterns of the stellar halo 
iRenda et allboOsHFont et al"200CT , an d the an gular momentum of disks (e.g. ICovernato et al.ll2006h . For 
a review of these issues, see .Brook et al.l (|2006al) . 



The ratio of a element to iron abundances has been used as a "clock" to infer different formation time- 
scales for different galaxy types, as well as for components of the Milky Way. Type II supernovae (SNe 
II), which trace star formation closely due the short life-spans of their high- mass progenitor stars, produce 
large quantities of oxygen, magnesium, silicon — the so-called a elements. By contrast, iron is produced 
predominately in Type la supernovae (SNe la), whose explosions are delayed. Hence, the high a-to-iron 
ratio characteristic of elliptical galaxies is interpreted as implying a rapid formation, while the low a-to-iron 
ratio of disks and dwarf galaxies implies a protracted formation. In our Galaxy, the low metallicity stellar 
halo, with peak iron abundance [Fe/H] —1.5 has a ratio [a/Fe] ~ 0.4, while the relatively metal rich 
thick disk (peak [Fe/ H] ^ —0.6) also ha s enhanced a elements compared w ith solar abundances. Recently, 
several studies (e.g. iBensbv et al.l 120051 : iReddv. Lambert. Sz Allendd 12006! ) have highlighted the difference 
in the a-to-iron element abundance ratio between the hot and cold disk compon ents of the Milky W ay, 
providing a clue to the processes involved in the formation of the components. In iBrook et al.l ( 2005[ l. we 
suggested that high star formation rates in gas-rich mergers result in increased metallicities with high a 
element abundances. Such events may be central in forming the thick disk component of the Milky Way. 
Here, we examine in detail a single isolated example of such a merger, and follow the [a/Fe] ratio, and how 
it evolves as the metallicity of the final galaxy increases during the merger event. 

Our study shows that two disk components naturally form in a gas-rich merger, and we examine struc- 
tural and kinematic features of the resulting galaxy. We show that the starburst which accompanies this 
merger has clear chemical signatures. We analyze and compare abundances in two specific groups of stars: 
merger stars, which are formed before and during the merger, and disk stars, those forming in the relatively 
quiescent period after the merger. The starburst which accompanies the merger results in rapid chemical 
enrichment with the large number of SNe II forming, ensuring high a element abundances. We discuss im- 
plications for the end products of gas-rich merger events, and the importance of such events in the formation 
of disk galaxies. This merger simulation will form part of a larger suite of simulations which explore the 
kinematic, structural, and chemical properties of gas-rich mergers remnants with various values of progenitor 
mass ratio, rotational velocities, and impact angles. 



Simulation Details and Results 



We simulate the merger using GCD-I-, which self-consistently models the effects of gravity, gas dynamics, 
radiative cooling, a nd sta r formation. We give an outline of the code here, whilst full details are found in 
Kawata fc GibsonI |2003h . GCD+ is a Tree/SPH algorithm that includes SNe la and SNe II feedback, 



and traces the lifetimes of individual stars, enabling us to monitor the chemical enrichment history of our 
simulated galaxies. Star formation occurs in a convergent gas velocity field where gas density is greater 
than a critical density, pcrit = 2 x 10~^^g/cm . The star formation rate (SFR) of eligible gas particles is 
then dp^/dt = dpg/dt — c^^pg/tg where c* — 0.5 is a dimensionless star formation efficiency, and tg is the 
dynamical time. This formula corresponds to the Schmidt law: SFR oc p^-^. The mass, energy, and heavy 
elements are smoothed over the neighboring gas particles using the SPH smoothing kernel. We assume that 
lO^^ergs is fed back as thermal energy from each SNe. Gas within the SPH smoothing kernel of SNe II 
explosions is prevented from cooling, creating an adiabatic phase for gas heated by such SNe. This adiabatic 
phase is assumed to last for the lifetime of the lowest mass star that ends as a SN II, i.e., the lifetime of an 
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8M0 star (100 Myr). This is similar t o a model presented in Ixhacker fc CouchmanI |2000l ). For SNe II, the 
metallicity-dependent stellar yields of IWooslev Sz Weaver ( 1995) are adopte d. For low- and intermediate- 
mass stars, we use the stellar yie l ds of Ivan den Hoek fc GroenwegenI (119971). W e adopt SNe la model of 



Kobavashi. Tsuiimoto. fc Nomotol ( 2000 ). and the yields of 



Iwamoto et al. 



U9M). 



The initial conditions are two galaxies with exponen tial gas disks, embedded in dark matter halos. These 
are created using GalactlCS (jKuiiken fc Dubinskilll995| ). and are essentially stable in the that their density 
profiles, potential, and velocity ellipsoids will not change significantly when individual galaxies are evolved. 
GalactlCS uses the lowered Evans model for the dark matter halo, which l eads to a constant-density core . 
Nunierical simulations suggest that halos have a cuspy central density profile ( Navarro. Frenk. fc Whitel[l996l 
1997t iMoore et allll999l : iGhigna et aDboOOl : Ijing fc Sutolbood iKlvpin et al.l boOlL However, these results 
are in conflict with many observations, so we prefer to use a constant-density core. In any case, our results 
should not be sensitive to the details of the central density profile, since the core contains only a small 
fraction of the total mass. 

The larger galaxy has total mass of 5 x IO^^Mq. The mass ratio is 2:1, disk scale lengths are 4.5 and 3.1 
kpc, and each has baryon fraction of 17%. Each galaxy consists of 40,000 baryonic and 100,000 dark matter 
particles. The small galaxy approaches with its orbital angular momentum vector inclined 17° to that of the 
large galaxy. Both galaxies have prograde rotation with respect to the angular momentum of the system, 
the orbital energy of the system is 1.7x10"''' ergs. The system is bound, and the spin parameter (ratio of 
orbital energy to binding energy) is A = 0.04. The gas disks evolve and form stars prior to the merger, and 



have gas fractions of fg i 
metallicity of log(Z/Zo) 



-^gas/(Afgas + A/stars) — 0.91 at the time of the merger. Gas is given an initial 
—4, and [a/Fe] = 0.35. We follow the simulation for 1.5 Gyrs. 



The gas-rich merger results in a final galaxy with disk morphology. Figure [T] shows the B-band luminosity 
m ap of the resultant ga l axy af ter 1.5 Gyrs of the simulation, where we employ the simple stellar population 
of lKodama. fc Arimotd (|l997l ). Shown, both face on (upper panels) and edge on (lower panels) are all stars 
(right panel), what we will call merger stars, which are those formed before and during the merger (middle 
panel), and disk stars, which are those that form after the merger (left panel). From the face-on view, 
it is evident that the remnant of this particular merger simulation is a ring galaxy, which indicates that 
prograde gas - rich d isk-disk mergers can produce ring galaxies. We find this especially interesting in light of 
Lavery et al.l (j2004l ). who found that the incidence of rings increases rapidly with rcdshift. We will study 
this further in our future papers. 

The B-band luminosity profiles are shown in Figure [2l and used to determine scale- lengths of 5.1 and 
4.1 kpc respectively for the merger and disk stars. A bulge component, within the inner ^ 3kpc, is apparent 
from the surface density profile, which is best fit with a bulge, disk, and thick disk. This is consistent 
with the merger remnant depicted in Figure 2 of Robertson et al. (200a). We note that our idealized initial 
disks do not have spheroid stellar components. Any stars in such component would end up in a spheroidal 
component in the final galaxy, and thus the idealized initial conditions are partially responsible for the lack 
of a significant spheroidal component in the final galaxy. Our results remain valid so long as any initial 
spheroid component is of low enough mass not to effect the dynamics of the merger. 

The star formation rate. Figure [31 shows a starburst which peaks at around 38OM0yr~^, during the 
merger. The end of this starburst is used to divide stars into merger and disk stars. Prior to the merger, star 
formation is around 30 MQyr"'^. After the merger event, the star formation rate drops below lOMQyr"'^ after 
around a Gyr, although this is affected by our idealized initial conditions that assumes a dense gas disk. The 
mass of stars born before, during, and after the merger, are 6.3 x 10^, 33 x 10^, and 19 x lO^M© respectively. 
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Fig. 1. — B-band Luminosity for final distribution, shown face on (upper panels) and edge-on (lower panels) 
of all stars (right), merger stars, (stars formed before and during the merger, middle), and disk stars, (formed 
after the merger, left). X, Y, Z are cartesian coordinates, with Z being along the axis of rotation. 

The merger stars have a thicker morphology than disk stars. In the disk region, 4 < Rxy < 10 kpc and 
\Z\ < Ikpc, the merger and disk stars have total stellar mass of S.SxIO^Mq, and 4.4 x 1O^M0 respectively. 
These two components also separate when we consider their metallicity distribution function. Figure [4] shows 
the metallicity distribution function for merger and disk stars. The merger stars have a peak metallicity 
of around [Fe/H] —0.8, while the disk stars peak at [Fe/H] ^ —0.2. The long low-metallicity tail of the 
merger component is due to a difference in metallicity between merger stars born before and during the 
merger- induced star burst. 

The top panel of Figure [5] shows the rotation velocity of merger and disk stars. Not surprisingly, the 
merger stars have significantly slower rotation. It is worth mentioning that a part of merger stars are counter- 
rotating. We also measure the rotation excluding the counter-rotating stars, and found that their rotation is 
still much lower than the disk stars. The timing of our division between merger and disk stars is consistent 
with the change in velocity dispersion, as seen in the bottom panel of Figure [5] which plots dispersion of 
rotational velocity versus formation time of the stars. 

The chemical abundances of the merger simulation are explored in Figure [6l The average abundances of 
O, Mg, and Si are used to calculate [a/Fe], and we take a slice through the galaxy at \z\ < Ikpc. The [a/Fe] 
versus radius plot (top left panel) shows the merger stars (x symbols) have [a/Fe] ~ 0.35 independent of 
radius (although the inner regions, associated with a bulge, is a little lower), whilst the later- forming disk 
stars (-1- symbols) have a value ^ 0.15 dex lower. The plot of [a/Fe] versus [Fe/H] in the top middle panel 
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Fig. 2. — B-band Luminosity profiles of merger stars (x symbols) and disk stars (+ symbols). Dotted lines 
are fits between 2.5 and 15 kpc, and indicate scale-lengths of 5.1 and 4.1 kpc respectively (ignoring a dense 
region associated with a ring in the disk stars). 
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Fig. 3. — Star formation rate (Mo/yr) against time in Gyrs. The starburst during the merger has a peak 
SFR of 380 M0 /yr. The dashed line shows the division between what we have termed merger stars and disk 
stars. The dotted line indicates where the merger begins. 
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Fig. 4. — The metallicity distribution function for stars in the region 4 < Rxy < 10 and \Z\ < Ikpc in the 
final galaxy. All stars are plotted as a solid line, and broken up into disk stars (dotted line), and merger 
stars (dashed line). 

is interesting, and motivated this letter. The merger stars maintain a higher ratio, [a/Fe], even when their 
[Fe/H] has increased to solar values, whilst the disk stars which are formed with values of [Fe/H] ~ —0.5 do 
so with relatively low values of [a/Fe]. Examination of the bottom left panel, where we plot the evolution of 
[a/Fe] with time, helps to explain this result. The first stars have [a/Fe] ^ 0.35, close to the initial condition, 
and their ratio drops prior to the merger. The value of [a/Fe] increases during the time of the starburst, 
as indicated by the jump from the second to the third points, before decreasing to the value found in the 
disk. During the starburst, associated large numbers of SNe II provide a supply of a elements which allows 
the maintenance of a high ratio of [a/Fe], even as the iron content of the gas reservoir is increased (bottom 
middle panel). After the starburst, the pollution of SN la becomes effective, and along with the contribution 
of abundances of low- and intermediate- mass stars, the value of [a/Fe] is maintained at a constant value. 
Further information on disk galaxy formation is found in vertical abundance gradients. The upper right 
panel of Figure [6] shows that neither merger nor disk stars have vertical gradient of [a /Fe] , while disk stars 
show a slight decrease in [Fe/H] at larger vertical height. 



Discussion 



In gas-rich mergers, two disk components naturally emerge. A thick disk consists of stars born before 
and during the merger, although w e note that stars born before t he merger can also end up in the stellar halo, 
which they did in the simulation of ISpringel fc Hernquisd (j2005f ) . A thin disk forms rapidly at the end of the 
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Fig. 5. — Top Panel: Rotational velocity plotted against radius for merger stars (x symbols) and disk stars 
(+ symbols). When a counter- rotating component is ignored, the remaining merger stars are shown as * 
symbols. Bottom Panel: Velocity dispersion in the direction of rotation versus formation time of final stars. 
The dashed line corresponds to the end of the merger, from Figure [S] 
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Fig. 6. — In all plots, merger and disk stars use x and + symbols respectively. Upper panels, [a/Fe] versus 
radius (left), [Fe/H] {middle), and vertical height (for stars with 4 < Rxy < 10, right). Lower left panel, 
[a/Fe] versus time. Lower middle panel, [Fe/H] versus time. Lower right panel, [Fe/H] versus vertical height 
(for stars with 4 < Rxy < 10). 



merger. The high star formation triggered in such events can also help explain the chemical abundances of 
the stars. The starburst results in rapid increase in metallicity (Fig. [6l bottom middle panel), with the short 
star formation timescale ensuring that stars formed are enriched largely by SNe II, and hence merger stars 
maintain a higher [a/Fe], even when their [Fe/H] has increased to solar values. The stars which form after 
this starburst epoch are formed with low values of [a/Fe], even those formed with low [Fe/H] ^ —0.5, and 
this quiescent epoch means that the stars form in a thin disk, with low velocity dispersion (Fig. [5l bottom 
panel). 

In brook et"ZI |2004bl) . we found that t he thick disk scale-length of o ur simulated galaxy was shorter 
than that of the thin disk. The observations of Yoachim fc Dalcanton ( 20061 ) find that thick disk scale-lengths 
are systematically larger than those of thin disks. In this merger simulation, the hot merger star population 
has an exponential profile with scale-length larger than that of the later-forming disk star population. This 
may favor a significant gas-rich merger being a feature of disk galaxy formation. Further simulations will 
determine whether old, hot disks with scale-length larger than those of young cold disks result for a wide 
range of merger parameters. Caution is required in interpreting the result, as subsequent "inside-out" thin 
disk growth through g as infall (not present in t he cur rent simulation) may increase the thin disk scale-length 



([Brook et al.l l2006b( ) . lYoachim &: DalcantonI (j2006l ) further find that low- mass disk galaxies have larger 



thick:thin disk mass ratios. They interpret this as evidence for the formation of the thick disk by direct 
accretion of stars, as progenitors of lower mass galaxies will more easily expel their gas from their small 
potential well prior to the merger. Yet observations suggest th at low mass galaxies are in f act more gas-rich, 
both at low redshift (e.g. ISchombert. McGaugh. fc Edeiil2001r ) and high redshift (|Erb et al. 2006), favoring 
gas-rich mergers as an interpretation of the higher thick:thin disk mass ratio in low mass gala xies. Also , 
the growth of the thin disk in small galaxies is perhaps regulated by their low densities (jPalcantonl 120061 ) 
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and supernovae feedback, i.e. the high thick to thin disk ratio of low-mass galaxies may be the result of less 
growth of a thin disk. 

Our idealized study does not include cold accretion from the intergalactic medium and infall processes, 
which are important in the birth of disk galaxies. Further, in the hierarchical structure formation scenario, a 
range of merger histories exists, and to explain the predominance of thick disk components in disk galaxies, 
one cannot rely on a single merger event. Rather, a few or several significant gas-rich merger events are 



likely to occur at an early epoch in the formation of a disk galaxy (IBrook et al.l 120051 ). A cold mode of 
accretion from filamentary structures also occurs in this scenario. But our study supports the view that 
the violent accretion of gas-rich galaxies is central in producing thick disk properties. By ignoring cold 
accretion, our simplified study highlights the effect of the starburst associated with such mergers on the 
chemical abundances of the forming stars, in particular the high a element abundances at high metallicities, 
and vertical abundance gradients. It is yet to be shown that cold accretion alone, perhaps resulting in a 
thin disk which is heated by an infalling satellite, or the dispersion of large star clusters ( Kroupa 2002t 



Elmegreen &: Elmegreen 2006[) . can reproduce such chemical signatures. Gas- rich mergers as the dominant 



process in forming thick disks provide a natural explanation of observed abundance patterns and gradients 
in the Milky Way disk components. A high-redshift epoch of gas-rich mergers is emerging as an important 
phase in the very birth of disk galaxies. 

The simulation was performed at the Laboratoire d'Astrophysique Numerique, Universite Laval. CB, 
SR, & HM are funded by the Canada Research Chair program and NSERC. DK is a JSPS Fellow. 
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